IEKP-KA/96-07 
KA-TP- 18-96 
|hep-ph/9607286 
July 7th, 1996 



Global Fits of the SM and MSSM 
to Electroweak Precision Data 

W. de Boer^ * A. Dabelstein^t W. HoUik^l 
W. Moslems U. Schwickerath^? 

1) Inst, fiir Experimentelle Kernphysik, Univ. of Karlsruhe, 
Postfach 6980, D-76128 Karlsruhe, Germany 
2) Inst, fiir Theoretische Physik, Univ. of Karlsruhe, 
Postfach 6980, D-76128 Karlsruhe, Germany 
3) Inst, fiir Theoretische Physik, Univ. of Miinchen, 
James Franck StrafJe, D-85747 Garching, Germany 



Abstract 

The Minimal supersymmetric extension of the Standard Model (MSSM) with light 
stops, charginos or pseudoscalar Higgs bosons has been suggested as an explanation of 
the too high value of the branching ratio of the boson into b quarks {Rf, anomaly). 
A program including all radiative corrections to the MSSM at the same level as the 
radiative corrections to the SM has been developed and used to perform global fits to 
all electroweak data from LEP, SLC and the Tevatron. The probability of the global 
fit improves from 8% in the SM to 18% in the MSSM. Including the 6 — 57 rate, as 
measured by CLEO, reduces the probability from 18% to 15%. In the constrained 
MSSM requiring unification and electroweak symmetry breaking no improvement of 
Rb is possible. 
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1. Introduction 



Present LEP data show a too high value of Rb {3Aa) and a too low value of Rc 
(2(t) where Rb(c) is the ratio Rb{c) = ^zo^bb{cc)/^zo^qq- In the past it has been shown 
by several groups that it is possible to improve Rb using supersymmetric models 
with light charginos, stops or Higgses, which yield positive contributions to the Zbb 
vertex [Q-IP]. In this paper we perform an equivalent analysis of all electroweak data 
both in the Standard Model (SM) and its supersymmetric extension (MSSM). The 
conclusion is that only a moderate increase in the probability can be found for the 
MSSM as compared to the SM, if all present limits on supersymmetric particles and 
the ratio 6 ^ 57 are taken into account. The analysis was performed using all actual 
electroweak data from Tevatron [|ll|, 0, LEP and SLC [0, the measurement 



of ^I^I^I^I^^^p^ from CLEO [T^ and limits on the masses of supersymmetric particles 



iHilT^ll, [T|,|0|. 

For low values of tan/3 the diagrams with charginos and right handed stops in 
the vertex loop are dominant, while for high tan/3 the exchange of the pseudoscalar 
Higgs between the outgoing 6-quarks, which is proportional to rrib ■ tan (3 becomes 
important too. It should be noted that light stops and charginos will contribute to 
the 6 — > 57 rate with an opposite sign as the SM diagram of the top - W boson loop, 
so the predicted 6 — > 57 rate can become easily too small, if the stops and charginos 
are light. Therefore, getting Rb and b ^ sj right is not easy, since Rb wants sparticle 
masses near the experimental limits, while b —>■ s'j needs sparticles well above the 
experimental limits or a sufficiently large tan (3. 

2. boson on-resonance observables 

At the Z boson resonance two classes of precision observables are available: 

a) inclusive quantities: 

• the partial leptonic and hadronic decay width Fjj, 

• the total decay width Tz, 

• the hadronic peak cross section ah, 

• the ratio of the hadronic to the electronic decay width of the Z boson: Rh, 

• the ratio of the partial decay width for Z ^ cc {bb) to the hadronic width, 

Rc, Rb- 

b) asymmetries and the corresponding mixing angles: 

• the forward-backward asymmetries Ap^, 

• the left-right asymmetries v4{^. 



• the r polarization Pr, 

• the effective weak mixing angles sin^^fj^. 

Together with the quantity Ar in the correlation of the W mass to the electroweak 
input parameters G^, Mz and aem? this set of precision observables is convenient for 
a numerical analysis of the supersymmetric parameter space. In the following the 
observables defined above are expressed with the help of effective couplings. 



2.1. The effective Z-f-f couplings 

The coupling of the Z boson to fermions / can be expressed by effective vector 
and axial vector coupling constants v{fp a{jj in terms of the NC vertex: 

Jnc = 1^ Uff - (^iffl^) . (1) 



where the convention is introduced : c^r = cos"^ 9w = 1 — = M^/M| Input 
parameters are the /i decay constant = 1.166392 x 10~^ GeV~^, uem = 1/137.036 
and the mass of the Z° boson Mz = 91.1884 GeV. The mass of the W boson is related 
to these input parameters through: 

71 ~ 2slyM^ ■ 1 - ArMssM ("em, Mw, Mz, rriu ...) ' 
where the complete MSSM one-loop contributions are parameterized by the quantity 
ArM5SAf|23- Leading higher order Standard Model corrections Q to the quantity 
Ar are included in the calculation. 

The effective couplings "wfyj, a{ff can be written as: 

viff = ^z{v^ + Avf + ZMQf) 

a{jf = ^z {af + AaO . (3) 
and are the tree-level vector and axial vector couplings: 

= ll-2Qfsl, , af = 4. (4) 
Zz, Zm are given in eq. (|l3). The complete MSSM one-loop contributions of the non- 



universal finite vector and axial vector couplings Av^ , Aa^ have been calculated p5[ 



together with the leading two- loop Standard Model contributions [p^, ^ . They are 



derived in the 't Hooft-Feynman gauge and in the on-shell renormalization scheme]^ . 
Fig. 1 shows the MSSM one-loop Z f f vertex correction diagrams. 
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Figure 1: MSSM one-loop Z ^ ff vertex correction diagrams. i,j,k — 1,..,2(4) 
are chargino, neutralino and sfermions indices. No particle permutations are shown. 




Figure 3: Z boson wave function renormalization. 
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The non-universal contributions can be written in the following way: 



Avf = F^'^' + AFv, 
Aaf = F^^' + AFa. 

The Standard Model form factors Fy^ corresponding to the diagrams of figs. 1 and 
2 can be found e.g in refs. ||2^, The diagrams with a virtual photon are listed for 
completeness in the figures. They are not part of the effective weak couplings but are 
treated separately in the QED corrections, together with real photon bremsstrahlung. 
The non-standard contributions are summarized by 

AFv = Y.^v +Vf^Zl + af6Z{ 

i 

AFa = Y.Fa +Vf6Z{ + af6Zl 

i 

where the sum extends over the diagrams of fig. 1 with internal charged and neutral 
Higgs bosons, charginos, neutralinos and scalar fermions, each diagram contributing 

to the Z f f -vertex. The self-energy diagrams of fig. 2 with internal neutral Higgs, 
chargino, neutralino and sfermion lines determine the field renormalization constants 

6Zl = -^v{m})-2m}[J:'y{m}) + E',{m})] 

SZ{ = S(m5) (5) 

with the scalar functions Tjv,a,s in the decomposition of the fermion self-energy ac- 
cording to 



The contributions from the Higgs sector are given explicitly in ref. [2S]. For the 



genuine SUSY diagrams, the couplings for charginos, neutralinos and sfermions are 



taken from |30], together with the diagonalization matrices given in section 3. 



The universal propagator corrections from the finite Z boson wave function renor- 
malization Zz and the 7Z mixing Zm are derived from the (7, Z) propagator matrix. 
The inverse matrix is: 

where S^, S^, S^^ are the renormalized self energies and mixing. They are obtained 
by summing the loop diagrams, shown symbolically in fig.|5, and the counter terms 
and can be found in ref. 
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The entries in the (7, Z) propagator matrix: 

/ A 



A 



fii/ 



7 ^"iZ 

A^z Az 



(7) 



are given by: 



fc2-M|+Sz(fc2) 



k^-M^ + ±z{k^ 



The renormahzation condition to define the mass of the Z boson is given by the pole 
of the propagator matrix (eq. P). The pole fc^ = M| is the solution of the equation: 



Eq. 



7^e [ (M| + S,(M|) ) Sz(k2) - E^z(M|) ] = . 
yields the wave function renormahzation Zz and mixing Zm'- 



(9) 



ResMz^z 



1 



l + 7^eS'z(k2) -7^e 



S2„ 
7Z 



(k2) 



k2+E^{k2) 



^,z{Ml) 



M| + S,(M|) 



(10) 



Z boson observables 
The fermionic Z boson partial decay widths F^j can be written as follows: 



ff 



where 



"qcd 



NcV2GM 



127r 

^^^^^ 



1 - Atmssm) 



+ 1.405(^)2 - 12.8(^ 



02 
4 7r2 



+ («/ 



/ = leptons 
/ = quarks 



(11) 



(12) 
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2) / 



bb 



NcV2G,Ml 



127r 



(13) 



In ATjji the b quark specific finite mass terms with QCD corrections p6[ are 
included. (Jq^'d given in eq. (^). 

The total decay width Tz is the sum of the contributions from leptons and quarks: 



/ 



(14) 



In the following Thad = J2q ^qq is the hadronic decay width of the Z boson. 
The hadronic peak cross section is defined as 

127r TeeThad 



Ml r| 



The ratio of the hadronic to the electronic decay width is defined as 

had 

J^p = — . 



(15) 



(16) 



The ratio of the partial decay width for Z ^ bb (cc) to the total hadronic decay width 
is given by 



bb(cc) 
had 



(17) 



The following quantities and observables depend on the ratio of the vector to axial 
vector coupling. The effective flavour dependent weak mixing angle can be written 
as 



sm 



The left-right asymmetries are given by 



1 + {viff/aij^Y 



while the forward-backward asymmetries can be written as 



Af^B = \A^A'. (20) 



3. The MSSM 



3.1. Higgs sector 

The scalar sector of the MSSM is completely determined by the value of tan /3 = 
V2/V1 and the pseudoscalar mass Ma, together with the radiative corrections. The 
latter ones are taken into account in terms of the effective potential approximation 
with the leading terms ~ m^, including the mixing in the scalar top system [^|. In 



this way, the coupling constants of the various Higgs particles to gauge bosons and 
fermions can be taken over from substituting only the scalar mixing angle a by 
the improved effective mixing angle which is obtained from the diagonalization of the 
scalar mass matrix. 

3.2. Sfermion sector 

The physical masses of squarks and sleptons are described by a 2 x 2 mass matrix: 

2_/ M^ + m} + Ml{li -Qfsl,) cos 2(3 mj(Aj + /i{cot /5, tan/?}) 
f ~ I mf{Af + fi{cotp,tanp}) M^^ + + MfQ/s^ cos2/? 

(21) 

with SUSY soft breaking parameters Mg, M^, Mj^, Af, and fi. It is convenient to 
use the following notation for the off-diagonal entries in eq. (pi]): 



A'j = Af + fi{cot (3, tan (3} . (22) 
Scalar neutrinos appear only as left-handed mass eigenstates. Up and down type 



sfermions in (^) are distinguished by setting f=u,d and the {u, d} entries in the 
parenthesis. Since the non-diagonal terms are proportional to mj, it seems natural 
to assume unmixed sfermions for the lepton and quark case except for the scalar top 
sector. The t mass matrix is diagonalized by a rotation matrix with a mixing angle 
^mix- Instead of Mq, M^j, Mfj, A[ for the b, t system the physical squark masses 
mi^,mi^, can be used together with A[ or, alternatively, the stop mixing angle 

For simplicity we assume m^^ = mi^, and tt, rf, c, s to have masses masses 
equal to the h squark mass. 

A possible mass splitting between hi-ti yields a contribution to the p-parameter 
p = 1 + Ap in terms of |]22|: H 



Apt, ^ + < - 2^^^ log ^ I . (23) 



"The superscript Ap" indicates that no left-right mixing is present. 
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As a universal loop contribution, it enters the quantity 

Ar~ AaEM-4^Ap + ... (24) 

and all the Z boson widths 

Fjj ~ 1 + Ap + • • • 

and is thus significantly constrained by the data on Mw and the leptonic widths. 
3.3. Chargino/Neutralino sector 

The chargino (neutralino) masses and the mixing angles in the gaugino couplings are 
calculated from soft breaking parameters Mi, M2 and /i in the chargino (neutralino) 
mass matrix[Q. The validity of the GUT relation Mi = 5/3 tan^ 9w M2 is assumed. 
The chargino 2x2 mass matrix is given by 

/ M2 MwV2smP \ 

= [ M^V2cos(3 -p j ' 

with the SUSY soft breaking parameters fi and M2 in the diagonal matrix elements. 
The physical chargino mass states xf are the rotated wino and charged Higgsino 
states: 

xt = v,,^l 

Xi = U.j^J ■ «,j = l,2. (26) 

Vij and Uij are unitary chargino mixing matrices obtained from the diagonalization 
of the mass matrix eq. |2^: 



U*A^^±V-i = diag(m-±, m^±) . (27) 



The neutralino 4x4 mass matrix can be written as: 
/ 



Ml -Mzsin^v^cos/? M^ sin ^vf sin /? ^^ 

M2 Mz cos 6w cos (3 —Mz cos 6w sin /3 

Mz sin 6y/ cos [3 Mz cos Ow cos /3 /i 

V Mz sin Ow sin /5 —Mz cos 6*54/ sin /? /i / 

(28) 

where the diagonalization can be obtained by the unitary matrix Nif 



N*A<^oN-^ = diag(m^o) . 



(29) 



The elements f/jj, Vij, Nij of the diagonalization matrices enter the couplings of the 
charginos, neutralinos and sfermions to fermions and gauge bosons, as explicitly given 
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in ref. | 
ref. 13 



30]. Note that our sign convention on the parameter fi is opposite to that of 



4. Results 



4.1. Chargino Masses 

As mentioned before the low tan (3 scenario of the MSSM needs a light right handed 
stop and light higgsino-like chargino for a large value of Rt, whereas in the high tan /3 
scenario one needs in addition a light pseudoscalar Higgs A [|1], |^. A higgsino - like 
chargino can be obtained for a low value of the parameter fi in the mass matrix (eq. 

Figs. ^ and ^ show the dependence of the chargino masses on the parameter /i in 
a region of the parameter space which yields a good global x^- case of high tan/3, 
m^2 is almost symmetric around zero, whereas in case of low tan/? this dependence 
is more complicated, as can be seen from fig. ^. For M2 = 3 \ n \ the light chargino 
mass passes zero at = —40, so the following low tan/? plots were made for /i > —40 
and /i < —40 GeV. The asymmetric structure of fig. ^ is reflected in the contours of 
constant Rb in the m^.^ versus light scalar top m^- plane (see fig. ||). High values of 
Rb up to 0.2194 are possible (see figs. |^ and |^, although these special regions of the 
parameter space are already experimentally excluded by the lower limits on sparticle 
masses. 

Taking M2 =| | does not change these results very much, as can be seen from a 
comparison of the distributions in fig. |^ (M2 = 3 | /i |) and fig. ^ (M2 =| /i |). The 
small increase of the at chargino masses around 80 GeV in the left hand part of fig. ^ 
is due to neutralino threshold singularities, for which an additional contribution 
has been added, if the sum of two neutralino masses is close to the mass. The 
sharp increase of the function at low chargino masses is due to experimental limits 
on chargino, neutralino and stop masses from LEP 1.5 [jT6|, 0, 0. 



4-2. Optimization of Parameters 

An optimization of free parameters of the MSSM was performed by minimizing 
a function using MINUIT [^. Several contributions to the x^ were taken into 
account: 

• experimental limits on the masses of supersymmetric particles and neutralino 
production from LEP 1.5 and Tevatron |T0, O, p 



precision measurements of on resonance observables from LEP [|14[, taking error 
correlations into account 

the measurement of the branching ratio sR{b^cep) from CLEO |15| 
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experimental limits 


m + 

Xl,2 


> 65 GeV 


•ffl 
^3,4 


> 13 GeV 

> 35 GeV 

> 60 GeV 
< 2 MeV 




> 48 GeV 


mh,mH,mA,mH± 


> 50 GeV 



Table 1: Mass limits assumed for the optimized fits. 



The experimental limits included in the fit are summarized in table ^ [jT], |Ty, |13, 
0, The calculation of the total decay width of the Z boson into neutralinos is 
based on reference the calculation of the ratio 6 57 on reference ^3 . 



As already mentioned, Rb mainly depends on the stop mass and the light chargino 
mass for the low tan j3 and on the pseudoscalar Higgs mass and chargino mass for the 
high tan (3 scenario. In order to get a feeling for the size of the effects, we study the 
dependence of the function on these parameters with optimization of the remain- 
ing parameters. This has been done for both the low and the high tan/3 solutions. 
The best solutions will be presented in the next chapter. 



• Low tan (3 : 



Fig. 10 shows the change in the best obtainable in the chargino - stop plane. 
For each value of the lighter scalar top mi and m^± in a grid of lOx 10 points an 
optimization of rrit, and the stop mixing angle ^mix was performed, assuming 
M2 = 3 I /i I for a fixed value of tan/3 =1.6. In the next section this assumption 
will be dropped. Low sparticle masses yield a sharp increase in the Ax^ in 
fig. |10| because of the included mass limits. The minimum is obtained for 
chargino masses just above the experimental limit, although it increases only 
slowly with increasing sparticle masses. Rb increases significantly with decreas- 
ing values of the stop and chargino mass, as can be seen from fig. |Tl|. Much 
less significant is the improvement of Re- Within the plane of fig. |TI| it changes 
less than 0.0005 units. The increase of Rb must be compensated by a decrease 
of as (see fig. |12D in order to keep the total Z-width constant. The stop mixing 
angle ^mix, shown in fig. 13, is mainly determined by the CLEO measurement 
of 6 — ^ S7. The chargino contribution to 6 — > 57 is proportional to the Higgs 
mixing parameter /i, which changes its sign for mj^± ~ 60 GeV (see fig.ffl), so 



the b ^ S'j rate changes rapidly for these chargino masses, as shown in fig. 14 



The uncertainty in the predicted 6 — » 57 rate from the renormalization scale 
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has been taken into account and was varied between mb/2 and 2mf, [^]. The 
scale itself has been chosen to be rrih- 



• High tan (3: 



Similar fits can be performed for the high tan/3 scenario in the pseudo scalar 
Higgs iriA versus light chargino plane. As in the low tan (3 case M2 = 3 | | 
was assumed. In fig. 15 the resulting change in the is given for fixed tan (3 = 
50. For small chargino masses there is a sharp increase in the due to the 
corresponding mass limit, see above. The best values for Rh can be obtained 
for small values of tua and m^^i, see fig. ^ As in the low tan/3 case the 
enhancement of Rb must be compensated by a decrease of a^, see fig. [l^, and 
the improvement of Rc is small, less than 0.0006 within the given parameter 
plane. The mixing angle,shown in fig. |l^, is mainly determined by the 6 — 57 
rate, which can be fitted in the whole rriA - chargino plane, see fig. 0. 



4-3. Best Solutions 

Standard Model Fits: 



The input values from Mw , rrit, the electroweak mixing angle and the line 
shape observables have been summarized in table H The SM predictions were ob- 



tained from the ZFITTER package and all the error correlations were taken from 
141. The fits were made with a,, nit and niH as free parameters, which resulted in 



as = 0.1215 ±0.0036 
mt = IQr.Stri GeV 



mn = GeV. 



The quoted errors have been determined using MINOS [T^. Further details of the 
procedure are described elsewhere, see for example ^ 0. The x^/d-o.f of the 
SM fit is 23.2/15 which corresponds to a probability of 8%. Here, the main contribu- 



tions to the originate from Rj, (Ax^ = 10.7), sin^ ©iTf from SLD (Ax^ = 3.6) and 



eff 



Rc (Ax^ = 3.3). The correlation parameter between tjih and mt for the best fit is 



approximately 0.7; this strong correlation is shown in fig. [2^. One observes that the 
upper limit on the Higgs mass is obtained for ^ 175 GeV; however, the upper limit 
is sensitive to sin^ ©1// shown by the dashed contour in fig. ^ where the precise 
value of sin e[ff from SLD was excluded from the fit. The dependence of sin^ ©1// 
on the SM Higgs mass is approximately logarithmic (see fig. ^1]). The LEP data 
alone without SLD yields tjih = U^lllt^ GeV, SLD alone yields mn = 15lf GeV, 
as indicated by the squares in fig. pll The latter value is excluded by the lower 
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limit of 58.4 GeV from the combined LEP experiments 



In addition, the ALEPH 
^ The Ax^ 



Collaboration gives a recent limit of 63.9 GeV on the SM Higgs mass |T9 
dependence of the Higgs mass is shown in fig. 



for various conditions. 



MSSM Fits and Comparison with the SM: 

In order to obtain the best MSSM fits the assumption M2 = 3 | // | is dropped and 
M2 is treated as a free parameter. Rb increases with decreasing tan [3. The fit results 
for a tan 13 as low as one are given in the first column of table the corresponding 
predictions of electroweak observables in table |^. Note the high value of and the 
corresponding low value of a^. The resulting /d.o.f. is 15.1/11, corresponding to a 
probability of about 18%. Unfortunately, the corresponding b ^ S'j rate is about one 
order of magnitude too small in this region of parameter space. Larger rates can be 
obtained either by heavier sparticle masses or by larger values of tan/3. With 6 — > 57 
included in the fit and a free tan/3, the preferred value of tan/3 was either around 1.6 
or 50. 

The fit results for these values of tan (3 are given in table |^ too and the predicted 
values of all observables with their pulls have been summarized in table ^. Note that 
the MSSM prediction of the W-boson mass is always higher than the Standard Model 
one. 

For the best solutions the gluino mass was fixed to 1500 GeV, the stau mass to 
500 GeV and the sbottom mass to 1000 GeV in both the low and the high tan/3 
scenario, since they are less sensitive to the LEP observables and therefore cannot be 
fitted. Their influence was studied by fixing them to different values and repeating 
the fits again. First the low tan (3 scenario will be discussed. A variation of the gluino 
mass from 200 GeV up to 2000 GeV did not change the best obtainable x^? varying 
the stau mass in the same range changed the less than 0.2. Furthermore, the best 
obtainable changed less than 0.2 when varying the sbottom mass and pseudoscalar 
Higgs mass from 800 GeV to 2000 GeV. For values of these two parameters below 
800 GeV the increased significantly, mainly because the prediction of R^ became 
too small. 

Within the high tan /3 scenario no significant change of the global x^ was detected 
when the gluino mass was varied between 200 GeV and 2000 GeV, but the stau 
mass was somewhat more sensitive. A variation of this parameter between 300 GeV 
and 700 GeV changed the global x^ less than 0.2, but if it was chosen higher than 
1000 GeV the x^ increased up to 1.6 units, mainly because the prediction of Ar 
became worse. The sbottom mass was varied between 800 GeV and 2000 GeV. As 
in the low tan/3 case there was no dependence on this parameter if it was chosen 
heavy, but for low values the preferred top mass became too small. M2 was fixed at 
1500 GeV. A variation between 1000 GeV and 2000 GeV changed the best reachable 
X^ less than 0.2, for smaller values the global x^ increased up to one unit. 

To check the influence of the assumptions on Mi for the best fits, the GUT relation 
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Ml = 5/3 tan^ 6w M2 was dropped and the fit was repeated with a free Mi, but this 
did not improve the best obtainable significantly. 

A direct comparison to the Standard Model fits is given for all three fits in 
figs. p^ - p6| . The resulting Standard Model x^/d-o.f. = 23.2/15 corresponds to a 
probability of 8%, the MSSM fits correspond to probabilities of 15% (tan/3 = 1.6, 
X^/d.o.f. = 16.9/12) and 10% (tan/3=50, x^/d.o.f. = 18.4/12). In counting the d.o.f 
the insensitive (and fixed) parameters were ignored. The high tan /3 scenario is hardly 
more probable than the SM, while for low tan f3 the probability of 15% is in between 
the best MSSM fit without 6 ^ 57 and the SM fit. 



Another interesting point are the predictions for as{Mz)- Fig. [23 shows a compar- 
ison of different measurements of as{Mz) (as given in pO|, ^) with the fitted values 
given in this paper. The fitted MSSM as{Mz) is slightly smaller than the Standard 
Model value and in better agreement with measurements from deep inelastic lepton 
scattering (DIS) and the world average [pO|] . 

4.4- Discovery Potential at LEP II 

The particle spectrum for the best fits, as shown in table |^, suggests that some 
SUSY particles could be within reach of LEP II. If they are not found, LEP II will 



provide stringent SUSY mass limits [pTj pSj . In the following the consequences of in- 



creased mass limits on the fits are discussed for both the low and high tan [3 scenario. 
Chargino Searches: 

The in the region of the best low tan/3 fit increases slowly for increasing 



chargino masses, see fig. ^ Chargino masses above a possible LEP II limit of 95 GeV 
will increase the global of the fit by approximately 2 units, which corresponds to 
a probability of 9%, which is hardly better than the SM probability, so one cannot 
consider the MSSM as a better solution in that case. For the high tan/3 scenario the 
dependence of the global x^ the chargino mass is small, see fig.^, so increased 
limits will hardly change the probability of about 10% for the best possible fit in this 
scenario. 

Stop and Neutralino Searches: 

The best tan/3 = 1.6 fit has a light stop mass of about 48 GeV, and the lightest 
neutralino is about 20 GeV, so this solution is just in between the regions of the 
parameter space which are excluded by stop searches at LEPl and the Tevatron 
IBf: LEP I limits for light right handed stops are about 45 GeV, while DO limits 



exclude 52(70) GeV< m^-^ < 92(87) GeV for m^-^ = 20 (40) GeV. A fit with increased 
neutralino mass limits of 45 GeV yielded a best solution of x^/ d.o.f. = 17.8/12, 
corresponding to a probability of 12%, which is worse than the solution presented 
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above, but still better than the Standard Model one. A similar /d.o.f. was obtained 
if the light stop was required to be heavier than 90 GeV. The reason for this can be 
found in the flat distribution in fig. This figure is similar to fig. |10|, but here 
M2 was treated as a free input parameter. 

The stop mass in the high tan (3 solution is 53 GeV, and the lightest neutralino is 
quite heavy, above 70 GeV, so there is no conflict with stop searches. Furthermore, 
the high tan (3 is insensitive to the stop mass, so an inreased mass limit does not effect 
the fit. 



Higgs Searches: 



The mass of the light scalar Higgs is a sensitive function of the top mass. If the 
top mass is below 180 GeV, the Higgs mass for the low tan {3 fit (see tab. |^) should be 
observable at LEP II, especially if one includes the second order corrections, which 
will lower the Higgs mass by 10 - 15 GeV Within the high tan/3 scenario 



both neutral Higgs bosons are light and have practically the same mass. Increasing 
the Higgs limits up to 90 GeV increases the /dof to 23.7/12, corresponding to a 
probability of only 2%, which is much worse than the SM fit. The steep dependence 
of the x^ oil the pseudoscalar Higgs mass, which is mainly caused by a too small value 
of is shown in fig. 0. 



5. CMSSM and Ri, 

In ref. [^] fits to low energy data have been performed within the constrained 
MSSM (CMSSM). In this case unification of gauge and b-r Yukawa couplings is 
assumed. Reproducing the large mass splitting in the stop sector, as given in tab. |^, 
needs a very artificial fine tuning of the few free parameters of the CMSSM, especially 
for the trilinear and Yukawa coupling in the stop sector, which can drive the diagonal 
elements of the stop matrix, eq. ^ apart. Note that the off-diagonal elements of 
this matrix are too small to generate a large splitting, since the left-handed stop 
is considerably heavier than the top, implying that one of the diagonal elements is 
considerably larger than the off-diagonal elements. 

In addition, problems arise with electroweak symmetry breaking, since this re- 
quires /i > M2, while requires < M2 for a significant improvement. In conclu- 
sion, within the CMSSM neither the high value of Rb nor the low value of Rc can 
be explained, so in this scenario the solution must be sought in common systematic 
errors for all experiments, like PDB branching ratios in the charm - sector |^ causing 
the too high value of i?b and too low value of R^- 



6. Conclusions 



The MSSM provides a good description of all electroweak data. Rh values up to 



15 



0.2193 are possible, but the 6 — > 57 rate is too low in this case. With all mass bounds 
and the 6 ^ 57 rate included in the fit, the best x^/d-o.f in the MSSM is 16.9/12, 
which corresponds to a probability of 15% as compared to the SM probabihty of 
about 8%. This best MSSM fit with light stops and light charginos is obtained for 
tan/9 = 1.6. Another solution with tan/? = 50 and light Higgses has a probability of 
10%, which is not much of an improvement over the SM. 

The enhancement of i?5 is compensated by a decrease in as{Mz) from 0.1215 in 
the SM case to 0.116 in the MSSM. The latter is in somewhat better agreement with 
precise measurements from DIS at low energies {as{Mz) = 0.112 ± 0.005). 

The best solutions predict chargino, stop and pseudoscalar Higgs masses which 
may be detectable at LEP II. The high tan/3 solution, requiring light Higgses, can 
certainly be excluded at LEP II, if no Higgses are found. For the low tan P scenario 
it will be more difficult to exclude the SUSY explanation of the too high value of Rb, 
since even for chargino and stop masses above 95 GcV moderate improvements are 
still possible. On the other hand, it has to be pointed out that the large splittings 
within the stop sector, which yield the best fit, are very difficult to obtain in a natural 
way within a constrained MSSM model, since they require a very special finetuning 
of the trilinear and Yukawa couplings. 
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Symbol 


measurement 


best fit observables 






SM 


MSSM 


tan /3 and pull 




pull 


1.0 


pull 


1.6 


pull 


50 


pull 


1V± ^ v_J \~ V 


Ql 1884+ n 0022 


91.1882 


0.089 


91.1884 


0.000 


91.1884 


0.000 


91.1884 


0.000 


ryfGeVi 


2.4964 ± 0.0032 


2.4973 


-0.284 


2.4956 


0.241 


2.4958 


0.202 


2.4952 


0.363 


"-'ft L^-^'-'J 


41 49 ± 078 


41.452 


0.483 


/II /I /I 1 

41.441 


0.634 


/IT A AO 

41.448 


A r /I A 

0.540 


/II A OA 

41.429 


A 'I'I'I 

0.777 




20.789 ± 0.032 


20.774 


0.456 


on TOTO 

20.7878 


A AO'7 

0.037 


OA '70'70 

20.7878 


A AO'7 

0.037 


OA 'VOAO 

20.7898 


A AO /I 

-0.024 




0.0171 ± 0.0011 


0.0164 


0.675 


0.0165 


0.577 


A A 1 /I 

0.0164 


0.626 


A A 1 r 

0.0165 


A TAT 

0.505 


Rb 


0.2211 ± 0.0016 


o.zioy 


6. All 


A 1 AO 

0.2193 


1.128 


A 1 OA 

0.2180 


1 A 1 

1.915 


A 1 OA 

0.2180 


1.967 


Rc 


0.1598 ± 0.0069 


0.1723 


-1.805 


1703 


-1 522 


1706 

\J . A. I \J\J 


-1 559 


1706 

L/. i 1 \J\J 


-1 565 




0.1002 ± 0.0028 


0.1034 


-1.145 


0.1040 


-1.366 


0.1038 


-1.286 


0.1043 


-1.453 


^FB 


0.0759 ± 0.0051 


0.0740 


0.374 


0.0743 


0.319 


0.0741 


0.347 


0.0745 


0.281 




0.842 ± 0.052 


0.9336 


-1.761 


0.9361 


-1.810 


0.9356 


-1.800 


0.9361 


-1.809 




0.6180 ± 0.091 


0.6680 


-0.550 


0.6684 


-0.554 


0.6682 


-0.552 


0.6685 


-0.555 




0.1394 ± 0.0069 


0.1477 


-1.196 


0.1482 


-1.270 


0.1479 


-1.235 


0.1479 


-1.229 




0.1429 ± 0.0079 


0.1477 


-0.605 


0.1482 


-0.667 


0.1479 


-0.636 


0.1484 


-0.712 


^n2ei7;((g^^)) 


0.2320 ± 0.0010 


0.2314 


0.562 


0.23138 


0.624 


0.23141 


0.593 


0.23133 


0.669 


Afw/ [GeV] 
1 - M^/Ml 


80.33 ± 0.15 


80.370 


-0.265 


80.417 


-0.583 


80.422 


-0.616 


80.452 


-0.814 


0.2257 ± 0.0047 


0.2232 


0.531 


0.2223 


0.727 


0.2222 


0.747 


0.2216 


0.869 


mt [GeV] 


175 ± 9. 


167.3 


0.858 


172.8 


0.239 


172.1 


0.32 


168.0 


0.776 


sin2ei7;(A^R)(SLD) 


0.23049± 0.0005 


0.23144 


-1.900 


0.23138 


-1.773 


0.23141 


-1.834 


0.23133 


-1.682 



Table 2: Measurements of the observables [14] and the predicted results of the fits with minimum x^- The pulls are 
defined by (measurement - predicted value) / error of the measurement . 



Fitted SUSY parameters and masses 




tan /?— 1 


tan (3- 


=1.6 


tan f3 


=50 




r~\r\ h ^ e o/ 

ILKJ U ' o J 


b ^ s'-f inc. 


■m J rip VI 

as 


173ih7 
0.1104±0.0043 


172±6 
0.1161±0.0038 


168±6 
0.1162±0.0039 


M2[GeV] 


25±8 




36±23 






/i[GeV] 


35±53 




42±9 




76±28 


m,-JGeV] 


48±5 




48±5 




53±40 


(fornix 


-0.163±0.115 


-0.203±0.091 


0.0021±0.0054 


mA[GeV] 










50±5 


Particle Spectrum 


m,-JGeV] 


^ 1 TeV 


m,-JGeV] 


48 


48 


53 


mg[GeV] 
m;-[GeV] 


1 TeV 
0.5 TeV 


m^±[GeV] 
m -l[GeV] 


91 




106 




1504 


81 




69 




76 


m^-o[GeV] 


15 




21 




73 


m^-o[GeV] 


35 




38 




79 


m^-o[GeV] 
m^-o[GeV] 


90 




97 




714 


102 




102 




1504 


mh[GeV] 


97 


110 




50 


mH[GeV] 


^ 1.5 TeV 






98 


mA[GeV] 


1.5 TeV 






50 


mH± [GeV] 


^ 1.5 TeV 






143 


M^[GeV] 


80.4174 




80.4224 




80.4520 


BRib^s-y) /ir|-4 
BR{b^ceu) 1 ^ 


(0.19) 




2.05 




2.30 


xVd.o.f. 


(15.1/11) 


16.9/12 


18.4/12 


Probability 


18% 


15% 


10% 



Table 3: Values of the fitted parameters (upper part) and corresponding mass spec- 
trum (lower part). The errors on the parameters are parabolic ones. The parameters 
given in the first column gave a minimum /d.o.f. of 15.1/11 for the LEP - observ- 
ables. Here b ^ S'j was not included in the fit, but the resulting b ^ sj rate is about 
one order of magnitude too small. Including b ^ s'y rate as measured by the CLEO 
Collaboration [|T5[ requires a higher value of tan/3, which reduces the best x^/d-o.f. 
to 16.9/12 (second column). On the right hand side the results of the optimization for 
tan /3 = 50 are given. The dashes indicate irrelevant parameters which were chosen 
high. 
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Figure 4: Dependence of the chargino masses on the parameter /i for M2 =\ fi \ 
and M2 = 3 I /i I for tan/3 = 1.6, ~ 0.117 and rhf^ ^ 60 GeV. No optimization 
of parameters was performed here. The shaded regions indicate chargino masses less 
than 65 GeV which are excluded by LEP 1.5 and chargino masses less than 96 GeV. 
It can be observed that for positive values of /i two light charginos are easier to obtain, 
if n and M2 have similar values. 
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Figure 5: Dependence of the chargino masses on the parameter fi for tan f3 = 50. In 
this case M2=1500 GeV was used and no optimization of parameters was performed. 
The shaded regions indicate chargino masses less than 65 GeV which are excluded 
by LEP 1.5 and chargino masses less than 96 GeV. The light chargino mass depends 
on yU, whereas the heavy chargino mass is dominated by the value of M2 and keeps 
approximately constant at 1500 GeV. For a fixed given chargino mass there are two 
possible solutions corresponding to /i < and /i > 0, respectively. 




Figure 6: Rb in the light stop versus hght chargino plane with M2 = 3 | /i | and 
tan/5 = 1.6. The upper part shows the solution with fi < —40 GeV, in the lower 
part the one with fi > —40 GeV is displayed. In the latter solution quite high values 
for Rb are possible, as can be seen in the figure. The dashed line in the upper plot 
indicates the 2a lower limit of Rb. 
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Figure 7: Rb in the versus light chargino plane with M2 = 3 | yU | for the high 
tan (3 solution, fi was chosen positive here. In this case choosing the opposite sign for 
fi doesn't change Rb- 

20 



Figure 8: Dependence of the absolute for A* > ~40 GeV (left side) and /i < 
—40 GeV (right side), using M2 = 3 | |. At the different points of the grid no 
optimisation of parameters was performed, but they were fixed to values near a good 
minimum. 
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Figure 10: The Ax^ in the hght stop and hght chargino plane for tan/? = 1.6. At 
each point of the grid an optimization of nit,as and the stop mixing angle (pmix was 
performed with /x > —40 and M2 = 3 | |, including the ratio b ^ sj. 
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Figure 11: Rb in the light stop and hght chargino plane. Optimization as in fig. |T^. 
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Figure 12: as in the light stop and hght chargino plane. Optimization as in fig. |T^. 




Figure 13: Stop mixing angle (pmix in the light stop and light chargino plane. It is 
mainly determined by the branching ratio b s'-f. Optimization as in fig. [1^. 
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Figure 14: 6 — 57 in the light stop and hght chargino plane. For chargino masses 
higher than 60 GeV (and /i > 0) the predicted value is close to the CLEO measure- 
ment of 2.32 ± 0.67 X 10~^. Optimization as in fig. llO. 
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Figure 15: The Ax^ in the pseudo scalar Higgs and hght chargino plane for tan/5 ~ 
50. For each given and light chargino mass an optimization of mt ^ ols rhii and 
the stop mixing angle ^mix was performed, including the b ^ rate and with the 
irrelevant parameter M2 set to 1500 GeV. 
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Figure 16: Rb in the pseudo scalar Higgs and light chargino plane for tan/? = 50. 
Optimization as in fig. 
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Figure 17: in the pseudo scalar Higgs and light chargino plane for tan/5 = 50. 
Optimization as in fig. 
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Figure 19: 6 — > 57 in the pseudo scalar Higgs and light chargino plane for tan (3 = 50. 
The prediction is close to the CLEO measurement of 2.32 ± 0.67 x 10~^ within the 
whole parameter space. Optimization as in fig. 
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Figure 20: A^;^ = 1 and Ax^ = 4 contour lines for all electroweak data including 
sin^ ©1// from SLD (continous line) and without it (dashed line). The stars indicate 
the best fits. 
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Figure 21: Dependence of the SM sin^ ®!;// the Higgs mass. The top mass 
rrit = 175 ± 9 GeV was varied within its error, as shown by the dashed band labelled 
SM (upper (lower) boundary mt=166(184) GeV). The SLD and the LEP measurement 
of sin^ ©1// cire also shown as horizontal bands. Fits to the electroweak data prefer 
nit ~ 170 GeV and light Higgs masses, as indicated by the squares for the separate 
LEP and SLD measurements, while the star is the result of the combined fit to SLD 
and LEP data. Clearly, the SLD value yields a Higgs mass less than the recents limits 
of 63.9 GeV by direct Higgs searches at LEP (shaded area) [p!9| . 
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Figure 22: Dependence of the SM Ax^ on the Higgs mass for a free top mass, taking 
all data (continous line), all data without the SLD measurement of sin^ ®e// (dashed 
line) and all data without Rb (dotted line). 
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Figure 23: The Ax^ in the region of the best fit in the hght stop and hght chargino 
plane for tan/3 =1.6. Here the constraint on M2 was dropped. At each point of the 
grid an optimization of rrit, M2, cts and the stop mixing angle (f)mix was performed with 
II > —40, including the ratio b ^ and the requirement I" z^neutraiinos < 2 MeV. 
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Figure 24: Resulting observables for the fit given in table ^ for tan/5 = 1.0. A 
significant improvement of Rb can be observed here. The ratio b ^ s'y was not 
included in this fit here. 
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Figure 25: Resulting observables for the fit given in table ^ for tan/3 = 1.6. was 
fixed to 1000 GeV, and the gluino mass were fixed to 1500 GeV. Including 6 — > 57 
in the fit it is still possible to improve the prediction of Rf, with Supersymmetry a 
bit. 
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Data / MSSM(tanB=50) 
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Figure 26: Resulting observables for the fit given in table ^ for tan f3 = 50. rhi, was 
fixed to 1000 GeV, M2 and the gluino mass were fixed to 1500 GeV. It is possible to 
improve the prediction of Rb with Supersymmetry even for high values of tan (3, but 
the result is not as good as for low values. 
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Figure 27: Comparison of different measurements of with the fit results. The 
data has been taken from and . 
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